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ABSTRACT 

The sole instrument on the upcoming Ice, Cloud, and Land Elevation Satellite ( ICESat-2 ) altimetry mission is 
a micropulse lidar that measures the time of flight of individual photons from laser pulses transmitted at 
532 nm. Prior to launch, the Multiple Altimeter Beam Experimental Lidar (MABEL) serves as an airborne 
implementation for testing and development. This paper provides a first examination of MABEL data ac- 
quired on two flights over sea ice in April 2012: one north of the Arctic coast of Greenland and the other in the 
east Greenland Sea. The phenomenology of photon distributions in the sea ice returns is investigated. An 
approach to locate the surface and estimate its elevation in the distributions is described, and its achievable 
precision is assessed. Retrieved surface elevations over relatively flat leads in the ice cover suggest that 
precisions of several centimeters are attainable. Restricting the width of the elevation window used in the 
surface analysis can mitigate potential biases in the elevation estimates due to subsurface returns at 532 nm. 
Comparisons of nearly coincident elevation profiles from MABEL with those acquired by an analog lidar 
show good agreement. Discrimination of ice and open water, a crucial step in the determination of sea ice 
freeboard and the estimation of ice thickness, is facilitated by contrasts in the observed signal-background 
photon statistics. Future flight paths will sample a broader range of seasonal ice conditions for further 
evaluation of the year-round profiling capabilities and limitations of the MABEL instrument. 


1. Introduction 

NASA’s second-generation Ice, Cloud, and Land 
Elevation Satellite {ICESat-2 ), currently planned for 
launch in late 2016, will provide observations to quantify 
the changes in ice sheets and sea ice, and key insights 
into their behavior (Abdalati et al. 2010). To achieve 
this, precise laser measurements of surface elevation, 
building on the capabilities of its predecessor ( ICESat-1 ), 
will be acquired to assess ice sheet mass balance and 
processes, as well as the time -varying thickness and 
volume of sea ice in the Arctic and Southern Oceans. 
ICESat-2 will also measure sea surface height in the ice- 
free sub-Arctic seas and provide large-scale vegetation 


Corresponding author address: Ron Kwok, Jet Propulsion Lab- 
oratory, California Institute of Technology, 4800 Oak Grove Dr., 
Pasadena, CA 91109. 

E-mail: ron.kwok@jpl.nasa.gov 

DOI: 10.1175/JTECH-D-13-00120.1 


biomass estimates through the measurement of vegeta- 
tion canopy height. Combining data from the ICESat-2 
mission with existing and forthcoming altimetry datasets 
will yield a 15+ year record of elevation change. 

Instead of the analog lidar used for ICESat-1 ( Abshire 
et al. 2005; Schutz et al. 2005), ICESat-2 will employ 
a photon-counting approach to obtain better measure- 
ment sensitivity with lower resource (power) demands 
on the satellite platform. A high repetition rate (10 kHz), 
low pulse energy laser at 532 nm in conjunction with 
sensitive detectors will measure the range of individual 
photons scattered from the surface. Six across- track 
beams will profile the ice surface, and for ice sheets 
the multiple beams address the need for unambiguous 
separation of ice sheet slope from elevation changes 
(Zwally et al. 2011). 

Photon-counting (PC), or micropulse, lidars have 
been employed for many years in atmospheric (e.g., 
cloud and aerosol) sounding applications (Spinhirne 
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1993; McGill et al. 2002), but the use of this type of lidar 
for high-precision surface altimetry is relatively new 
(e.g., Harding et al. 2011). Whereas traditional altime- 
ters have used analog detection with fast digitization of 
the return waveform, PC lidars provide the distance 
traveled by individual photons. Surface elevations are 
estimated from the observed photon distributions, 
sometimes referred to as “photon clouds.” PC imple- 
mentations reduce the thermal and power demands for 
a satellite-based instrument, but a small field of view is 
required to limit the confounding effects of solar back- 
ground in the return data. For ICESat-2 , the diameter of 
each laser footprint on the surface is —10 m (compared 
to the ICESat footprint of —70 m) with a receiver field of 
view of — 40-m diameter. The design uses six transmit 
beams in three pairs in a rectangular configuration. A 
yaw angle from the satellite track then sets the specific 
spacing of the footprints (http://icesat.gsfc.nasa.gov/ 
icesat2). 

The Multiple Altimeter Beam Experimental lidar 
(MABEL) is an airborne PC lidar used as a technology 
demonstrator for the photon-counting instrument [Ad- 
vanced Topographic Laser Altimeter System (ATLAS)] 
on the ICESat-2 mission (McGill et al. 2013). Data from 
this instrument will be used to validate instrument 
models and to assess performance of ATLAS that will 
be launched on ICESat-2. Table 1 compares the per- 
formance of ICESat-2 and MABEL. The first Arctic 
deployment of MABEL in April 2012 acquired sea ice 
data with two dedicated flight paths, one covering 
>3000 km in length (8 April) and the other >1000 km 
(10 April) (Fig. 1). These survey paths were selected 
to provide a broad sampling of ice conditions prior 
to the expected onset of melt over the Arctic sea ice 
cover. The first ER-2 flight on 8 April (operations 
flown out of Keflavik, Iceland) acquired data over 
the Lincoln Sea — the compact, deformed ice north of 
the coast of Greenland — and the sea ice in the Fram 
Strait and east Greenland Sea between Greenland 
and Iceland. The second flight resurveyed the track 
just east of the Greenland coast (see Fig. 1). On 
10 April, the ER-2 was joined by the P-3 from Oper- 
ation IceBridge (Koenig et al. 2010). The P-3 carried 
the Airborne Topographic Mapper (ATM) instru- 
ment (Krabill et al. 2002), an analog lidar, to provide 
near-coincident coverage for assessment of MABEL’s 
performance. 

Since this is the first acquisition of this type of PC data 
over Arctic sea ice, the objective of this paper is to in- 
vestigate the phenomenology and to examine the utility 
and achievable precision of the MABEL instrument for 
retrieval of ice and open water elevations for freeboard 
calculations (e.g., in Kwok et al. 2009). The intent is to 


Table 1. Comparison of MABEL and ICESat-2 instrument 
parameters and expected-observed performance. 


Parameter 

ICESat-2 

MABEL-Channel 6 
April 2012 

Operational altitude, R (km) 

490 

20 

Wavelength, A (nm) 

532 

532 

Telescope diameter, At 

0.8 m 

12.7 cm 

PRF 

10 kHz 

5 kHz 

Pulse energy, E t (/jlJ) 

41/160 a 

1.4/0.2 b 

Footprint ( He 2 ) (mrad) (m) 

31 (15) 

100 (2) 

FOV (mrad) (m) 

83 (41) 

210 (4.2) 

Filter width (p.m.) 

30 

-150 

System-detector efficiency, 

3 

5 

ToQe (%) 

Swath width (km) 

±3.00 

±1.05 

Signal levels 0 over winter sea 

Expected 

Expected b Observed d 

ice (photons per shot) 

Snow-covered ice 

1.6/6.2 a 

1.4/0. 19 b 0.2 d 

(albedo = 0.9) 

Open lead (albedo = 0.15) 

0.26/1. 0 a 

0.2/0.03 b 0.04 d 


a Weak/strong beam. 
b Before/after fiber damage. 

c Lidar equation used to compute expected signal levels: 
N r = oLx(E t \lhc)(TAATToQElTrR 2 ), where a\ = surface albedo, 
T a = molecular atmospheric transmission (0.81 is used in 
calculations), Qe = detector quantum efficiency, T a = system 
optical transmission, c = speed of light, and h = Planck’s 
constant. 

d Observed signal levels after fiber damage (see sections 3 and 4). 

provide an examination of MABEL acquisitions over 
sea ice but not a comprehensive analysis of the PC in- 
strument for surface profiling. The paper is organized as 
follows. Section 2 describes briefly the MABEL in- 
strument, its data products, and the ATM acquisitions 
used in our analysis. In section 3, we examine the phe- 
nomenology of photon returns from the sea ice cover 
with segments of MABEL data acquired during this 
Arctic deployment. An approach to locate the surface 
and estimate its elevation in the photon clouds is de- 
scribed in section 4. The achievable precision and quality 
are assessed over several large leads. In section 5, we 
compare the retrieved elevations from MABEL with 
two near-coincident profiles from the ATM lidar. Sum- 
mary remarks and conclusions are provided in the last 
section. 

2 Data description 

In this section, we provide brief descriptions of the 
two lidar systems (MABEL and ATM) and their ele- 
vation data products used in this paper. MABEL data- 
sets and documentation are available online (http:// 
icesat.gsfc.nasa.gov/icesat2/data/mabel/mabel_docs. 
php). Datasets from Operation IceBridge (OIB) are 
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Fig. 1. MABEL flight paths (on 8 and 10 Apr) plotted on an Environmental Satellite ( Envisat ) synthetic aperture 
radar (SAR) mosaic of the region acquired on 8 Apr. Near-coincident ATM data from OIB were acquired on 10 Apr. 
Inset shows the ER-2 aircraft. MABEL is located in the nose of the platform. The segments discussed in the paper are 
from flight paths south of the Fram Strait. 


archived at the National Snow and Ice Data Center 
(NSIDC). In particular, the dataset of interest in this 
paper is the surface elevation data acquired by the Air- 
borne Topographic Mapper (ATM), an analog lidar 
used to provide high-precision lidar mapping of surface 
elevation. 

a. MABEL — Photon-counting lidar 

For this first deployment to the Arctic, the MABEL 
instrument was mounted in the nose of the ER-2 air- 
craft. The relevant instrument parameters are shown in 
Table 1. The MABEL laser has a pulse width of —2 ns, 
and the pulse repetition frequency (PRF) can be varied 
between 5 and 25 kHz. With a PRF of 5 kHz (used 
during these flights) and a nominal ER-2 ground speed 
of —200 ms -1 , a pulse is transmitted approximately 
every 4 cm along track. From an operating altitude of 
20 km, the telescope field of view (FOV: — 210/rrad) 
covers a — 4-m-diameter footprint on the ground and the 


laser illuminates a smaller spot of —2m (1/e 2 ) in di- 
ameter; this ensures that the illuminated spot is within 
the FOV of the telescope. An inertial measurement unit 
(IMU) is mounted directly on the telescope assembly to 
permit accurate determination of instrument pointing. 
A digital camera system (DCS) with a resolution of 
— lm captures aerial images of the surface along the 
MABEL tracks. These images are invaluable for iden- 
tifying the surface types (open water, snow-covered ice, 
thin ice, etc.) sampled by the lidar. 

The output of the MABEL laser at 1064 and 532 nm is 
split into eight 1064-nm and sixteen 532-nm beams. 
Returns from the two wavelengths facilitate the exami- 
nation of elevation biases due to subsurface scattering 
from water, and ice and snow volumes. In MABEL, each 
beam path can accommodate a small neutral density 
filter for varying the energy levels at each footprint. 
Coupling a given beam to one of 107 fibers (at each 
wavelength), which are pointed at a fixed set of look 
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angles, allows selectable off-nadir sampling to maximum 
view angles of ±3° (±52mrad). From a flight altitude of 
20 km, the beams are capable of acquiring data to 
±1.05 km on each side of the track (see Fig. 2 in McGill 
et al. 2013). With a separate set of transmitter fibers for 
each wavelength, the same surface area can be illumi- 
nated with both wavelengths (i.e., the footprints for each 
wavelength can overlap). For the two flights over sea ice 
of interest here, MABEL was configured with a narrow 
measurement swath. The 532-nm channels were set to 
illuminate the following look angles (relative to nadir) at 
approximately 0, ±0.1, ±0.2, ±1.0, ±1.9, ±2.0, ±3.0, 
±4.0, and ±5.0mrad. The 1064-nm channels were set at 
±0.1, ±0.2, ±1.9, and ±2.0 mrad. At the 20-km mea- 
surement altitude, this provided a cross-track width of 
±100m. 

The sixteen 532-nm receiver fibers are routed to 
a 16-channel Hamamatsu model H7260 photomultiplier 
tube (PMT) detector and the eight 1064-nm receiver fibers 
are routed to eight individual Excelitas single -photon- 
counting modules (SPCMs). For these detectors, there is 
a dead time ( t d ) of —2.5 ns. Detector dead time is time 
required for the detector to recover from a photon event 
(the triggering of the detector) before the next photon 
can be detected. For MABEL, the surface return and 
background rates are much lower than 400 MHz 
thus, an insignificant fraction of photons events occur 
during the dead time, so the sampled rate is linear with 
the input rate. Detected photons are time tagged for 
range determination and have a measured resolution of 
83 ps (about 12 mm). 

Analyzed MABEL products (release 7) contain geo- 
located elevations of individual photons between -1 
and 5 km relative to the World Geodetic System 84 
(WGS84) ellipsoid, a commonly used datum for mea- 
surement of elevations. Products are distributed as data 
files, each containing 60 s of surface returns from 300 000 
pulses, covering —13 km along track. Henceforth, we 
refer to these as MABEL file segments. 

b. MABEL instrument performance 

An unexpected physical deterioration of optical fibers 
in the transmit path over the course of this deployment 
affected each beam differently. The deterioration was 
due to damage sustained during the alignment of the 
laser and beam splitter prior to the beginning of the 
flights. This resulted in an overall degradation in instru- 
ment performance and data quality because of reductions 
in transmitted energy and received signal strength (i.e., 
reduced photon count from the surface). Prior to fiber 
damage, the MABEL signal levels were expected to be 
1.4 and 0.2 photons per shot (see Table 1) over a high- 
reflectance snow-covered sea ice surface (albedo: —0.9) 


and water in open leads (albedo: —0.15), respectively; 
the goal was to match the per-shot signal performance of 
the ICESat-2 instrument. 

To assess the impact of the fiber damage on signal 
levels, we selected one beam (beam 6) at the center of 
the array with the best signal strength. Based on our 
estimates, the fiber damage caused at least a 7-times 
reduction in pulse energy (from 1.4 to 0.2 julJ) with 
a similar reduction in the received signal level. Ex- 
pected per-shot photon count decreased from 1.4 to 
0.19 for snow-covered sea ice and from 0.2 to 0.03 for 
open water leads (see Table 1). These estimates can be 
compared with the observed returns of 0.2 and 0.04 
photon per shot from these two surface types (see Table 1; 
also discussed in sections 3 and 4). The near corre- 
spondence between the expected and observed signal 
levels are probably fortuitous because there are un- 
certainties in the albedo of the surfaces that we selected 
for comparisons. 

Even though the signal levels are reduced, signal 
photons from consecutive MABEL shots can be com- 
bined to simulate the signal levels of ICE Sat-2. How- 
ever, there are two considerations if that approach were 
used: first photon bias and background counts. Since the 
average surface signal per shot from high reflectance 
surfaces in the current MABEL dataset is less than one, 
the first photon bias in surface elevation estimation 
(described in Yang et al. 2011) is zero. This is good in 
practice, but this does not allow us to examine the im- 
pact of first photon bias in sea ice altimetry. Second, 
when multiple MABEL shots are combined to simulate 
the signal from a single ICESat-2 shot, the effective 
background rate (per shot) increases linearly. There are 
—17 MABEL shots (pulse interval: —4 cm) in an ICESat-2 
pulse interval (70 cm). If 17 shots were combined, 
a MABEL background rate of 0.1 MHz (estimated in 
section 3) would give an effective rate of 1.7 MHz. This 
is, in fact, close to what one expects to see in ICESat-2 
data from the polar regions. Also of note is that the 
degraded MABEL signal density (i.e., counts per meter) 
of 1.0-5.0 photons m _1 over leads/snow-covered ice is 
roughly equivalent to the expected ICESat-2 signal 
density of 1. 4/8.6 photons m _1 over these same surfaces. 
In any case, these are important factors to consider when 
examining this particular dataset. 

c. ATM lidar 

For a more detail description of the ATM instrument, 
the reader is referred to Krabill et al. (2002). Briefly, the 
ATM is a conical-scanning laser ranging system oper- 
ated at a wavelength of 532 nm with a pulse repetition 
frequency of 5 kHz and a scan rate of 20 Hz; the off-nadir 
scan angle is 15°. To provide accurate elevation estimates, 
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the instrument includes a differential GPS system for 
accurate positioning of the aircraft, and inertial sensors 
(accelerometers and gyroscopes) for precise attitude 
determination. With nominal OIB flight parameters 
(i.e., operating altitude and ground speed: 500 m and 
250 kt; lkt = 0.51 ms -1 ), the ATM scanning geometry 
provides an across-track scan swath of —250 m. Near the 
center of the swath, the spacing between neighboring 
laser footprints, of — lm in diameter, is approximately 
3-4 m in the along- and across-track directions. The 
sample density is higher (submeter) near the edges of 
the swath due to the conical scanning geometry of the 
system. We note here that this lidar does not provide 
contiguous spots on the ground. 

Generally, the backscatter from a snow or ice surface 
within the beam of the ATM is sufficient for measure- 
ment of the time delay of a return signal and for de- 
termination of the total propagation distance. However, 
the presence of extremely smooth surfaces along the 
flight path results in some measurement dropouts, due 
to weak surface returns at off-nadir angles. In the pro- 
cessing of the data, the travel time data are combined 
with GPS navigation measurements and aircraft orienta- 
tion parameters to derive surface elevation measurements 
relative to the WGS84 reference ellipsoid. Typical ele- 
vation accuracy is better than 10 cm (Krabill et al. 2002). 
Retrieved elevations are provided in data files that cover 
tracks of —35 km in length, each containing over a mil- 
lion elevation estimates. 

3. Phenomenology 

To carry out high-precision surface profiling of the sea 
ice cover with the PC approach, a first step is to un- 
derstand the expected returns from various surface 
types and contributing noise sources (system and geo- 
physical). Even though the physical basis of the lidar 
returns from the surface are not different from that of an 
analog lidar, there are important distinctions and con- 
siderations in the use of these data. With several ex- 
amples in this section, we show the signatures of sea ice 
returns within an elevation window close to the surface. 
As mentioned above, all the MABEL data used here are 
from channel 6 (532 nm). Since the dataset was acquired 
in early April, the sea ice conditions can be considered 
to be that of early spring; that is, there is very little 
surface melt except perhaps near the ice edge. Prior to 
melt, a significant fraction of Arctic sea ice is snow 
covered (of high albedo), while the remaining areas are 
covered by newly formed thinner ice (of lower albedo) 
that have yet to acquire a snow cover and by open 
water. Significant contrasts between these surfaces are 
expected. 


a. Signal versus noise photons 

One distinction between PC and analog lidars is the 
discrete nature of the noise. The total number of pho- 
tons (NP tot ) within a range window includes those 
backscattered from the surface (NP surf ) and atmosphere 
(NPatmX those from solar background (NP S bkg), and 
those introduced by noise in the detector (NP det : dark 
counts), namely, 

NP tot = NP surf + (NP atm + NP Sbkg + NP det ) . 

Clearly, the larger the fraction of surface (or signal) 
photons compared to those photons (background) from 
sources within the parentheses (i.e., higher signal-to- 
noise ratio), the better for surface retrieval. Background 
or noise photons add to uncertainties in the surface 
profiling and retrieval processes. On a cloud-free day, 
NP atm can be considered to be small compared to the 
other terms. For MABEL, the detector noise rate (at 
0.02 kHz) is much smaller than the signal rates and 
solar background discussed below. Henceforth, NP atm 
and NP det will be assumed to be negligible and not 
discussed. 

The solar background rate ( B s , usually measured in 
photons per second or MHz) is the solar zenith radiance 
due to surface and atmospheric scattering of solar en- 
ergy at the laser wavelength (A) that is detected by the 
instrument. The solar zenith angle varies latitudinally, 
seasonally, and with time of day. The quantity B s is 
negligible at night but significant when the sun is high in 
the sky. Neglecting atmospheric effects, this quantity 
can be calculated as 

B s — S\ a \ cos(0 z )(<S) (MHz) . (1) 

In the equation, B s is dependent on the bidirectional 
reflectance of the surface (a A ), the solar zenith angle 
(0 Z ), and the solar flux at the top of the atmosphere (, S A ). 
The quantity S includes constants and system parame- 
ters that are not germane to the current discussion. This 
calculated background rate (in MHz) when divided 
by 0.5 X (speed of light) gives the expected number 
of background photons found within a 1-m elevation 
window. 

The observed and expected [using a A = 0.8 in Eq. (1)] 
solar background rates along the two flight paths are 
shown in Fig. 2. The observed background rate is cal- 
culated using the total number of photons between an 
elevation of 4 and 5000 m above the surface in each 
MABEL file segment (or 300000 shots). The observed 
background rates are less than expected, but for long 
stretches of the 3000-km flight line on 8 Apr north of the 
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Fig. 2. Calculated (black line) and observed (solid circles) solar background noise rates over flight paths on (left) 8 
and (right) 10 Apr. See Fig. 1 for location of these flight paths. 


Greenland coast, the calculated and observed background 
noise rates are highly correlated. There is little vari- 
ability in albedo in this region of relatively compact 
snow-covered sea ice, and the changes in the solar zenith 
angle explain a large fraction of the variability in the 
solar background rates. This can be compared to the more 
variable rates over the mixed ice cover in Nares Strait 
and those just south of Fram Strait. Here, the mix of ice 
types from old ice to that of recently opened leads, with 
large contrasts in albedos, contributes to the observed 
variability. Similar variability can be found in the shorter 
flight line on 10 April (a repeat survey of the track flown 
on 8 April) of the mixed ice cover (deformed with rel- 
atively wide leads) just south of Fram Strait. The dis- 
crepancy between the calculated and observed solar 
background is due to unmodeled changes in albedo and 
to residual biases in system parameters that are cur- 
rently being addressed by the MABEL instrument team. 

For the two flight paths, the magnitude of B s (—0.1 MHz) 
translates into an approximate mean noise density of 
—3.3 X 10 -4 photons m -1 . This is equivalent to 1.5 X 
10 -3 photons over a limited elevation window contain- 
ing the surface (say, 5 m). The likelihood of finding 
a noise photon is relatively small compared to the 
number of expected surface (or signal) photons of be- 
tween 0.1 and 1.0 photon per shot from a dark sea ice 
lead or a snow-covered ice surface (see discussion be- 
low). It should also be noted that the detector dead time 
is much smaller (by three orders of magnitude) than the 
expected interval between background photons ( 1/B S ); 
thus, the impact of the background photons on the sur- 
face elevation retrieval is small. 

More importantly for sea ice, since B s provides 
a measure of surface albedo, it is a useful parameter for 
supporting the classification of surface types in the 
photon data (see discussion in next section). 


b. Sea ice in two MABEL file segments 

We use two examples (see Figs. 3 and 4) to show the 
signatures of different surface types in the geolocated 
photon clouds from MABEL. Of particular geophysical 
interest are the contrasts in the signal and background 
rates for discrimination of ice and open water in the re- 
trieval of sea ice freeboard. Each example shows a 5. 5 -km 
subsegment extracted from a 13-km file segment. A 
large open water lead is centered in the first example 
(Fig. 3) while more along-track spatial structure in the 
ice cover is seen in the second example (Fig. 4). 

The DCS imagery (in Figs. 3a and 4a) shows the sea 
ice cover profiled by MABEL. Both examples contain at 
least four sea ice types as suggested by their image in- 
tensities and elevations (relative to the dark lead). Vi- 
sually, the brightest samples seem to belong to that of 
snow-covered sea ice; the next brightest is thin ice with 
a thin layer of snow; then thin ice that is bare of snow; 
and last, of open water in leads. We recognize the thin 
ice categories as thin ice because their elevations are 
closest to that of open water (see Figs. 3c-e and Figs. 4c-e). 
The derivation of surface elevation estimates shown in 
Figs. 3e and 4e are discussed in the next section. 

The red and black dots (in Figs. 3b and 4b) show the 
along-track surface (-3 to 4 m) and background noise (4 
to 5000 m) counts, respectively — henceforth, the fol- 
lowing notation is used to describe elevation intervals: 
[-3, 4] m and [4, 5000] m. In both examples, the profile 
of the background counts is positively correlated to the 
intensities of the DCS samples (p — 0.98). This is ex- 
pected as the PC system is similar to a camera system that 
records scattering of incident solar flux at 532 nm. The 
ratio of background (B) photon counts NP^ ow /NP^ ad 
from snow-covered ice surfaces and open water leads are 
of —5 — approximates what one expects from the albedo 
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Fig. 3. Returns from the sea ice cover (channel 6) along a cloud-free ~5.5-km MABEL track, 
(a) DCS image containing the MABEL track, (b) Total photon counts from two elevation 
intervals: surface (—3 to 4 m) and background (4-5000 m); photon counts are sums of returns 
from 125 shots (—5 m along track); note the different scales for the (left) background and 
(right) surface (right) counts, (c) Photon elevations; colors represents total counts for 125 shots, 
(d) Number of shots needed to aggregate at least 100 signal photons, (e) Retrieved elevation. 
Location of this MABEL segment is shown in Fig. 1. Transparent light blue band shows 
characteristics of photons originating near the sea surface. 


differences [of —0.8 for ice and — (0. 1-0.2) for leads] 
between the two surface types (Grenfell and Perovich 
1984; Perovich et al. 2002). 

In these examples, the surface or signal counts and 
the DCS image intensities are only correlated over the 
brighter surface categories of snow-covered sea ice, snow- 
covered thin ice, and thin ice. Over the dark sea ice leads 
with low background counts, it is striking that the signal 


or surface photon counts are nearly 5-6 times the surface 
counts from a snow-covered ice surface (see Figs. 3b 
and 4b). This is contrary to what one expects from a 
Lambertian surface, and it suggests that these high sig- 
nal counts are associated with near-nadir specular re- 
turns from very smooth surfaces in the open lead. This 
contrast in surface photon density can be seen in the 
color-coded plot of photon elevations in Figs. 3c and 4c, 
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Fig. 4. As in Fig. 3, but for a different segment. 


where the colors represent 125-shot photon counts in the 
neighborhood (centered) of each detected photon. How 
close to nadir do these returns originate? From consid- 
eration of the attitude and relative alignment of the 
IMU and the lidar, the returns seen in this and examples 
that follow are within l°-2° of nadir. 

Away from the specular returns, the ratios of surface 
(, S ) photon counts NP^ ow /NP^ inice between the snow- 
covered ice surface and thin ice in both the solar back- 
ground (B) are ~(3-4). As expected the contrast is 
somewhat lower than the albedo differences (of —0.8 for 
ice and >0.2 for thin ice) between snow-covered ice and 
open water. 


A measure of PC performance or the sensitivity of the 
system is the number of shots it takes to accumulate 
a certain number of near-surface or signal photons (an 
aggregate). Here, we use 100 as the photon count (which 
can be converted to shots per aggregate, where the ag- 
gregate count = 100). For MABEL over sea ice, the 
average varies between —100 and 1400 shots per ag- 
gregate (see Figs. 3d and 4d), or from 1 photon per shot 
(near specular surfaces) to 1 photon per 14 shots (low 
reflectance surfaces). This parameter is fairly stable over 
snow-covered sea ice (—500 shots per aggregate), but 
there is considerable variability in thin ice and open 
water areas (shaded regions in Figs. 3 and 4). Thus, even 
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though there are no surface returns from certain shots, 
the MABEL photon clouds provide a nice depiction of the 
surface profile of the sea ice surface because of the high 
pulse repetition rate or small sample spacing (— 4 cm). 

c. Elevation distributions 

Here, we examine the elevation distributions of three 
distinct surface types: lead, lead with specular return, 
and snow-covered sea ice (Fig. 5). The character of these 
distributions informs the design of surface elevation 
retrieval and classification procedures. 

In the first example, we show the elevation distribu- 
tion of a sea ice lead with low signal and background 
counts (Fig. 5, left panel). The surface and background 
counts (Fig. 5a), and the associated photon cloud of the 
— 1.45-km-wide lead (identified by a horizontal black 
line in Fig. 5b) are shown within the context of a 13-km 
MABEL file segment. An exploded view of the photon 
counts and elevation cloud are shown in Figs. 5c and 5d, 
and the elevation distribution of the surface photon 
cloud can be seen in Fig. 5e. From the 35 982 transmitted 
laser pulses over the —1.45-km-wide lead, there were 
only 1591 surface photons within the [— 2, 1] m elevation 
interval in the neighborhood of the surface. This gives 
approximately one signal per surface photon in —22 
shots. From the elevation plot (Fig. 5d, left panel), we 
see a higher density of the photons below the surface on 
the right edge of the lead but noticeably higher back- 
ground counts from the same location. The extended tail 
in the elevation distribution of the return from the ag- 
gregate of 1591 photons (see Fig. 5e, left panel) suggests 
returns from the subsurface [elevation standard de- 
viation (SD) = 0.40 m], but it could also be due to sur- 
face roughness. As noted below, the returns may also be 
skewed as a result of the laser pulse shape. 

The next example (see Fig. 5, middle panel) shows 
specular returns from a — 500-m-wide lead that is 
flanked by several lower reflectance leads. The signal- 
surface counts from the specular lead are almost 5 times 
higher. The background counts from this specular lead, 
which are comparatively lower than those from adjacent 
leads, suggest a smooth surface of ice or water. In this 
case, the signal and solar background counts do not 
covary because of the specular return from the surface. 
There are 7470 signal photons in the 12497 shots or — 1 
signal photon in 1.7 shots; this can be contrasted with the 
1 -in-22 ratio computed for the lead in the previous ex- 
ample. With the larger number of counts from a smooth 
specular surface over a short distance, the apparent 
elevation distribution from this photon aggregate (see 
Fig. 5) provides a measure of the shape of the transmit 
pulse. The width of the return (SD = 0.17 m) is nar- 
rower than that of the previous example and compares 


well with the expected pulse width of 2 ns (SD = 15 cm) 
of the laser, which has a sharper leading edge with 
a trailing edge that decays at a slower rate. In these two 
examples, we note that this contrast in both elevation 
and surface/background photon counts (as seen in this 
and the next example) is extremely useful for the dis- 
crimination of ice and open water in freeboard calcula- 
tions. 

The 2-km segment shown in Fig. 5 (right panel) is from 
a highly reflective snow-covered surface. Both the sur- 
face and background counts are consistently higher than 
those seen in the left and middle panels of Fig. 5. 
Compared to the photon elevations of the leads in the 
previous examples, their variability is also higher. In the 
1.6-km window, there is —1 signal photon in 5 shots, 
higher than the counts from leads but lower than that 
from a specular surface. The width of the elevation 
distribution from the photon aggregate is 0.54 (SD), 
characteristic of returns from a surface that is rougher 
than those discussed above. As noted earlier, since the 
average surface signal per shot from high reflectance 
surfaces in the current MABEL dataset is less than one, 
the first photon bias (described in Yang et al. 2011) is 
zero and does not need to be considered. 


4. Surface elevation retrieval 

In traditional altimetry, the elevation of a surface is 
estimated from a digitized waveform from the surface 
scatter. For a PC system, one approach is to create an 
elevation distribution (or “waveform”) from the photon 
clouds (like those seen in Figs. 5a and 5b) and then use 
traditional approaches to refine the surface elevation. 
One advantage of the high PRF and small spot size of 
PC systems is that there are a number of ways one could 
select photons to construct these waveforms for surface 
elevation estimation, which are adapted to the length 
scale of the surface of interest (e.g., the width of a sea ice 
lead). In this section, we first describe an approach to 
retrieve surface elevations in the photon cloud. Second, 
the key parameters used to control the quality of the 
retrieval process are discussed. Finally, examples are 
provided to illustrate the effectiveness of this process 
with MABEL data. 

a. An approach 

We assume that the expected return, s e (h), from a 
surface with a Gaussian elevation distribution (width w 
or 2 times the standard deviation) at an elevation offset 
of h 0 [i.e., G(h\h 0 , w)] can be written as 

s e (h\ h D , w) = s t (h) * G(h\ h Q , w) , 


( 2 ) 
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Fig. 5. Elevation distributions of (left) leads, (middle) open water, and (right) snow-covered surfaces (channel 6). (a) Photon counts 
along a cloud-free MABEL file segment (300000 shots and —13 km long; note the different scales): surface (-3 to 4 m) and background 
(4-5000 m) and 125-shot photon counts, (b) Photon elevations; colors represent 125-shot photon counts (—5 m along track), (c) Photon 
counts in the subsegments indicated by the dark black line in (b). (d) Photon elevations within the subsegments, (e) Elevation distributions 
within the subsegment (bin size = 2.5 cm); quantities in the top-right corners show the distribution means and standard deviations (in 
parentheses). The plots in the (a),(c) middle panels have the y-axis labels of the left and right panels. 
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where s t (h) is the impulse response of the system and * is 
the convolution operator; that is, we assume the ex- 
pected return to be the convolution of the system im- 
pulse response with a Gaussian distribution. 

Further, we define the mean squared difference be- 
tween the normalized received signal (or elevation dis- 
tribution) s YQC (h ) and the normalized expected return, 
s e (h), to be 

++ w ) = ix ~ ++ H ’)l 2 

W N 

X S 'rcc (A) = 1 

N 

'Zs e ( h i -, h 0 ,w) = l. (3) 

N 

The two variables s e (h) and s YQC (h) are normalized be- 
cause we have no a priori information about the am- 
plitude of the received signals (or the reflectivity of the 
surface). 

Then, we consider the minimum in the error surface, 
e 2 (h 0 , w), to provide the best estimate of the height 
offset, h Q , of a surface with a Gaussian distribution of 
width w. This can be represented as 

(h 0 ,w) = argmin e 2 (h Q ,w) h g € [h v h 2 ] w e [w v w 2 ]; 

( 4 ) 

that is, the location of the given argument for which the 
function e 2 attains its minimum value within the in- 
tervals defined by h 0 E [hi, /q] and w E [uq, uq] is 
considered the best estimates of h Q and w. 

b. Key parameters in retrieval process 

The key parameters that control behavior of the re- 
trieval process described above are the number of signal 
photons and bin size used in constructing s YQC (h); the 
window size or the number of bins [N in Eq. (3)] of the 
elevation distribution used in the process; and the size of 
the search space in width (w E [uq, uq]) and elevation 
offset ( h 0 E [hi, /q]). 

An elevation distribution s YQC (h) can be constructed 
using photon elevations from a fixed or variable number 
of shots. Rather than using photons from a fixed number 
of shots, we find it best to use a fixed number of photons 
(or aggregates) over a variable number of shots. Over 
surfaces with lower returns, waveforms contain photons 
from a larger number of shots, and vice versa. This al- 
lows the control of signal strength (or signal-to-noise 
ratio) of the constructed waveform and thus the con- 
sistency in the quality of the surface retrievals. This is 
well suited for sea ice surfaces with its broad range of 


expected reflectance. We note, however, that the spatial 
resolution of the retrievals is proportional to signal 
strength and varies along track. Below, we show exam- 
ples of retrieved elevations using different aggregate 
counts. As for bin size, we use a width of 2.5 cm, which 
is approximately the elevation resolution of the PC 
system. 

The window size specifies that a portion of the ele- 
vation distribution be used in the surface estimation 
process. We use the standard deviation (cr) of the dis- 
tribution within —2 to 3 m of the mean elevation to de- 
fine the window size. In the results shown here, we 
truncate the distribution at 2a from the mean; that is, we 
use a window size of 4cr. This effectively limits the 
number of noise photons used in the retrieval process. 
Potentially, narrowing the width of the window could 
also be used to reduce the sensitivity of the retrieved 
elevation to the contributions of subsurface returns 
present in the tails of the distribution. For the dataset 
examined, the [-2, 3]m interval works quite well; 
however, we may have to expand the window to include 
returns from taller ridges. 

To implement our approach, we also require esti- 
mates of s t (h), that is, the impulse response of the sys- 
tem. Surface retrievals are sensitive to the shape of s t (h), 
especially when the surface is relatively flat compared to 
the pulse width. Asymmetric transmitted pulse shapes or 
system response, if not accounted for, would introduce 
biases in the elevation estimates. As measurements of 
the system impulse response are not available at the 
time of this writing, we use the elevation distribution of 
the specular return (Fig. 5, middle panel) in the previous 
section as a replica of the impulse response pulse. This 
model of the replica, if not contaminated by subsurface 
return, suggests asymmetry in the impulse response. 

The search space in elevation offset and width seen in 
Eq. (4) (i.e., w E [uq, uq] and h Q E [hi, h 2 \) depends on 
the quality of the initial estimate of the mean elevation 
calculated above and the expected distribution of 
the surface relief. For the results discussed below, 
the following intervals are used: h Q E [-0.5, +0.5] m, 
w E [0.0, 1.5] m. The mechanization of the retrieval 
process is not discussed here, suffice it to say that 
finding the minimum in e 2 (h 0 , w) can be computa- 
tionally optimized. Efficiency can be obtained by 
pregenerating s e (h ) over the search space of interest. 

c. Retrieval examples 

Here, we discuss the surface retrieval results from 
a 6.5-km MABEL segment using aggregates of 100 
photons to construct the elevation distributions s rec . 
Figure 6 shows the photon cloud (Fig. 6a), the estimated 
half-width of the Gaussian (G) (Fig. 6b), and the 
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Fig. 6. Surface elevation retrieval using 100-photon aggregates, (a) Photon cloud from a 6.5-km segment of 
MABEL, (b) Half-width of Gaussian from retrieval process, (c) Retrieved elevation, (d) Elevation distribution and 
modeled return (red) [at locations A-C in (c)]. (e) Error surface with the light and dark gray lines showing the local 
minima along the x axis (half-width of Gaussian) and y axis (elevation offset) and with colors showing the relative size 
of the minima. Panels (a)-(c) have the same x axis as in (c). 
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retrieved elevation profile (Fig. 6c). Over a relatively flat 
portion of the ice cover near the center of the segment 
(with low photon density — bluish color), the half- width 
of the Gaussian varies between 0 and 0.3 m, and the el- 
evation varies by —8 cm. This can be compared to the 
more deformed (Gaussian half -width ranges up to 1.5 m), 
snow-covered surfaces (of higher photon density) that 
ranges up to 2 m in elevation above the lowest point in 
the profile. 

In particular, at three locations (A, B, and C in Fig. 6c) 
along the MABEL segment, we show the elevation 
distributions [s Y&c (h)], the modeled return, [, s e (h ) in red], 
and the error surface e 2 (h 0 , w) used in the estimation of 
the surface elevations (Figs. 6d and 6e). These locations 
were selected to illustrate the effectiveness of the re- 
trieval process over a range of roughness scales. The 
distribution in red can be thought of as the best fit of the 
modeled returns [s e (h)] to the s YQC (h). 

The s Yec (h) at A (Fig. 6d, left panel) is that of a smooth 
surface with a relatively narrow distribution that ex- 
hibits the asymmetry, with a noticeable tail, that is seen 
in the specular return in Fig. 5e (left panel), that is, our 
assumed replica of the system impulse response. The red 
line in the figure shows the modeled return fitted to the 
distribution. A very sharp (i.e., the local curvature; Fig. 6e, 
left panel) and well-defined minimum is evident in 
e 2 (h 0 , w). The s Yec (h) at B (Fig. 6d, middle panel) is that 
of a broader return. In this case, as the surface roughness 
broadens the return distribution, the expression of the 
asymmetry in s t (h ) is no longer apparent. And, the 
minimum in e 2 (h 0 , w) (Fig. 6e, middle panel) is not as 
sharp as that seen in A (Fig. 6e, left panel) but never- 
theless clearly identifiable in the search space. The 
s Yec (h) at C (Fig. 6d, right panel) is that of an even 
broader return with a bimodal character that suggests 
photons from a mixture of smooth and deformed ice 
(perhaps ridges of different heights). Similarly, the ex- 
pression of the asymmetry in s t (h) is no longer apparent. 
The minimum in e 2 (h 0 , w) (Fig. 6e, right panel) is 
broadened further than that seen in A or B. Still, the 
minimum is evident but the surface is less symmetric 
about the minimum: the curvature of the surface in the 
direction that varies the width of the Gaussian is re- 
duced. This is likely due to the mixture of surface 
roughness in the distribution, as the surface behavior 
departs from that of the model assumption. 

These examples illustrate the importance of the pulse 
shape in the retrieval of elevation over smooth surfaces. 
As the elevation distribution deviates from that of 
a Gaussian when a mixture of ice types are present, the 
minimum in e 2 (h 0 , w) becomes less pronounced. The 
curvatures may provide an indicator, in addition to 
the width of the Gaussian, of the character of the surface. 


d. Sensitivity to aggregate count 

As discussed above, the quality and sensitivity of re- 
trieved elevations are sensitive to the size of the photon 
population (aggregate) used to estimate the mean ele- 
vation. Two examples (in Figs. 7 and 8) from MABEL 
illustrate the differences in the retrieval results using 
aggregates of 50 and 100 photons. Both examples in- 
clude flat smooth surfaces and rougher snow-covered 
surfaces. The smoother surfaces are ideal for assessing 
the relative precision of the instrument and retrieval 
procedure, and the rough surfaces are useful for un- 
derstanding the impact of spatial sampling. 

The photon clouds in Figs. 7a and 8a show the distri- 
bution in elevation as well as the photon densities in 
elevation (in terms of counts) along the track. There are 
clear distinctions between lower surface counts from 
relatively smooth surfaces and higher counts from 
rougher, snow-covered surfaces. The contrasts between 
the elevation profiles using 50 and 100 photon aggre- 
gates can be seen in Figs. 7b and 7c, and 8b and 8c. Not 
unexpectedly, the retrieved elevations using 50-photon 
aggregates seem more sensitive to surface detail and 
thus noisier overall. For MABEL, the average spatial 
resolution of the elevation profile is ~8 and —16 m for 
the 50- and 100-photon aggregates, respectively; spatial 
resolution depends on the number of shots needed to 
accumulate the signal photons used in the aggregate (as 
indicated in Figs. 7d and 8d). 

The second moment (standard deviation) of the sur- 
face elevation profiles is useful in assessing the impact of 
aggregate sizes. In Figs. 7e and 8e, we show the standard 
deviation of the elevation profile over 100-m intervals. 
Over relatively flat areas in Fig. 7, we find standard de- 
viations of —(5-6) (black) and 2-3 cm (red) in the 50- 
and 100-count aggregates, respectively. This suggests 
that an along-track precision of several centimeters can 
be obtained with 100-photon aggregates with this par- 
ticular implementation of the MABEL. Another ex- 
pected result is that the standard deviation of elevations 
in the rougher ice is somewhat more muted using the 
100-count aggregates because of smoothing of the ele- 
vations when larger populations are used. In general, the 
retrieval process seems well behaved when the aggre- 
gate counts are varied. 

5. Comparisons with profiles from the ATM lidar 

On 10 April the MABEL (on ER-2) and ATM lidars (on 
the P-3 from Operation IceBridge) flew near-coincident 
ground tracks to obtain surface elevation profiles for 
assessment of the MABEL surface elevation retrievals. 
The ATM instrument on the P-3 with an approximate 
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Fig. 7. Example comparison of retrieved surface elevations using 50- and 100-photon aggregates: (a) 125- shot 
photon counts, (b) 50-photon aggregate elevation, (c) 100-photon aggregate elevation (d) shots and spacing per 
50- and 100-photon aggregates, and (e) 100-m elevation standard deviation for 50- and 100-photon aggregates. 


scan swath of —250 m, an analog lidar also at 532 nm, 
increases the likelihood of coincident coverage (Figs. 9a 
and 9b). 

The elevation profiles from the two instruments 
for two different — 4-km segments are shown in Fig. 9. 
These comparisons are produced in two steps: the cal- 
culation of average profiles of the two datasets and 
the location of the MABEL ground track in the ATM 


swath. Given the approximate resolution of the MABEL 
surface retrieval (— 15 m), we produce surface elevation 
estimates from ATM at approximately the same length 
scale. It should also be noted neither the ATM nor 
MABEL instruments provide contiguous ground sampling. 
In the case of ATM, there are gaps associated with 
the scanning processes and for MABEL there are spots 
with no surface returns. Elevations for comparison with 
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Fig. 8. As in Fig. 7, but for a second example. 


MABEL retrievals were calculated by computing 15-m 
(diameter) Gaussian averages of the ATM elevations 
centered on each MABEL location. 

Accurate collocation of MABEL tracks in the ATM 
swaths is challenging because of ice drift and the dif- 
ferent ground speeds of the instrument platforms (ER-2 
and P-3). Ice movements and time delays between ac- 
quisitions limit our ability to register finescale surface 
features along track. As the sea ice moves and deforms, 
time differences in coverage introduce misregistrations 
between the surface profiles, especially when movement 


of discrete ice floes are especially evident at the meter- 
level length scales considered here. In the two examples 
shown in Fig. 9, the delays between the two instruments 
are —18 min. The impact of ice motion is appreciable. At 
a speed of >10 cm s -1 (typical for this area), sea ice drifts 
—140m in 20 min but the strain rates (or relative motion 
between floes) are less predictable. Hence, the expected 
misregistration between surface features in the presence 
of surface motion, even with meter-level geolocation ac- 
curacy, is not small. The comparisons (in Fig. 9) are pro- 
duced by maximizing the correlations between the two 
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Fig. 9. Comparisons of surface elevation profiles from MABEL with near-coincident ATM profiles, (a) (top) ATM 
swath with MABEL ground track (black) and (bottom) ATM (red) and MABEL (black) profiles. The time dif- 
ference shown in the top left corner of the (bottom) panel is between the two airborne platforms, (b) As in (a), but for 
a different segment along the track. (c),(d) Scatterplot of elevation samples in ATM and MABEL profiles in (a),(b), 
respectively; numerical values of correlation and differences between the profiles are shown in the top-left corner of 
each panel; elevations are averages of ATM footprints within 15 m (diameter) of MABEL locations. 


elevation profiles (along and across track) around a best 
initial location of the MABEL ground track within the 
ATM swath. As there are residual biases in absolute ele- 
vation (relative to the WGS84 reference ellipsoid) in both 
instruments, the two profiles are tied together (i.e., relative 
biases removed) at the points indicated by the arrows (in 
Fig. 9) once the offset with maximum correlation is found. 

Even with considerations in spatial coregistration 
and sampling discussed above, the two comparisons 


(scatterplots in Fig. 9c) between MABEL and ATM 
surface profiles nevertheless give correlations of 0.9 
and 0.86 with similar standard differences of 0.3 and 
0.26 m, respectively. Broadly, the results show that the 
MABEL and ATM retrievals covary along track and 
thus are sensitive to similar surface features. The scatter 
is symmetric around the line with unity slope, which 
suggests that thickness-dependent biases between the 
two retrievals are not evident — at least in this dataset. 
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However, detailed analyses of the differences may not be 
not warranted given the contributions of the different 
error sources discussed above. 

6. Conclusions 

In this paper, surface elevation profiles of the Arctic 
Ocean sea ice cover acquired by the MABEL in- 
strument are examined. Over the next several years, 
MABEL will serve as a tool to support the development 
of algorithms for retrieval of geophysical parameters 
during the ICESat-2 mission. Since this is the first Arc- 
tic deployment of this type of instrument, we provide 
a broad survey of the phenomenology, one approach to 
retrieve surface elevation, and an assessment of the re- 
trieval approach. Specifically, this paper addresses the 
sea ice dataset acquired on 8 and 10 April 2012. Even 
though damage to the optical fibers during this campaign 
resulted in degraded performance of the instrument, the 
acquired dataset provided a valuable first look of the 
Arctic sea ice cover with a micropulse lidar. In this 
section, we summarize the results as well as outline fu- 
ture plans on the use of MABEL to provide a broader 
spectrum of seasonal ice conditions (winter and sum- 
mer melt) for further evaluation of the profiling capa- 
bilities and limitations of the MABEL and ATLAS 
instruments. 

For the ICE Sat-2 mission, the ATLAS instrument will 
be tasked to provide basinwide sea ice freeboard esti- 
mates for calculation of ice thickness. MABEL will 
support the demonstration of this requirement to carry 
out high-precision surface profiling of the sea ice cover 
with a photon counting (PC) approach. Although there 
are important distinctions and considerations in dealing 
with the photon clouds from PC systems, the physical 
basis of the lidar is not different from that of an analog 
lidar. We examined the surface signatures of sea ice 
returns, and the contrasts between the photon distri- 
bution and signal strength from different snow-covered 
sea ice types and various forms of thin ice and open 
water. Specular returns, from what appears to be thin 
ice or open water, are particularly distinctive in the 
MABEL returns. Subsurface or rough surface returns 
seem to be present in the photon clouds in some of 
the MABEL file segments over sea ice leads. But, the 
unambiguous identification and likelihood of such oc- 
currences are more difficult to quantify. This remains 
a topic of investigation with upcoming MABEL data 
acquisitions. 

An approach to retrieve surface elevations in the 
observed photon cloud was introduced. Elevation dis- 
tributions of the surface are constructed using a fixed 
number of photons (or aggregates) over a variable 


number of shots and then compared with modeled 
returns. Fixed-sized aggregates allow the control of 
signal strength and the consistency in the quality of the 
surface retrievals. Returns are modeled as the convo- 
lution of the system impulse response with that of 
a Gaussian distributed surface elevation distribution. 
This is well suited for sea ice surfaces with its broad 
range of expected reflectance. Restricting the width of 
the window used in the surface analysis can mitigate 
potential biases in the elevation estimates due to sub- 
surface returns at 532 nm. At this writing, the impact of 
subsurface return on elevation retrievals awaits the 
analysis of larger volumes of MABEL data from future 
acquisitions. Retrieved surface elevations over relative 
flat sea ice leads suggest that precisions of several cen- 
timeters seem attainable. 

On 10 April, the MABEL (on ER-2) and ATM lidars 
[on the P-3 from Operation IceBridge (Koenig et al. 
2010)] flew near-coincident ground tracks to obtain 
surface elevation profiles for assessment of surface ele- 
vation retrievals from MABEL. The two comparisons 
between MABEL and ATM surface profiles gave cor- 
relations of 0.9 and 0.86 with similar standard differ- 
ences of 0.3 and 0.26 m, respectively. The results show 
that the MABEL and ATM retrievals covary along- 
track and thus are sensitive to similar surface features. 
However, detailed assessment of the differences in this 
particular dataset may not be warranted given the con- 
tributions of the different error sources: resolution, 
sampling, registration, ice motion, and the variability of 
surface relief at the 10-15-m length scale. 

Over the next several years, MABEL wifi be tasked to 
provide data to support development of operational 
algorithms for retrieval of geophysical parameters. 
During the ICESat-2 mission, the ATLAS instrument 
wifi acquire year-round data over the ice-covered ocean 
of the Arctic and Antarctic, which includes all types of 
surface and atmospheric conditions during different 
seasons and through seasonal transitions. The sensitivity 
of surface retrievals to solar background noise and atmo- 
spheric effects such as signal reduction through absorption 
and scattering as well as multiple-scattering-induced range 
delay must be characterized. Understanding these geo- 
physical and system characteristics is crucial for the de- 
sign of robust retrieval algorithms. The goal, for MABEL, 
is to provide a broad spectrum of seasonal ice conditions 
(melt and summer) for further evaluation of the pro- 
filing capabilities and limitations of the MABEL and 
ATLAS instruments. Additionally, the issue of sub- 
surface returns at 532 nm on surface elevation retrieval 
remains. At this writing, the next deployment of MABEL 
wifi be to acquire data over the summer ice cover of the 
Arctic Ocean. 


1168 


JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY 


Volume 31 


Acknowledgments. The digital camera system (DCS) 
on the ER-2 is operated and maintained by the Airborne 
Sensor Facility under the auspices of NASA’s Airborne 
Science Program. RK and GFC carried out this work at 
the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aero- 
nautics and Space Administration. 

REFERENCES 

Abdalati, W., and Coauthors, 2010: The ICESat-2 laser altim- 
etry mission. Proc. IEEE, 98 , 735-751, doi:10.1109/ 
JPROC. 2009. 2034765. 

Abshire, J. B., X. Sun, H. Riris, J. M. Sirota, J. F. McGarry, S. Palm, 
D. Yi, and P. Liiva, 2005: Geoscience Laser Altimeter System 
(GLAS) on the ICESat mission: On-orbit measurement 
performance. Geophys. Res. Lett., 32 , L21S02, doi:10.1029/ 
2005GL024028. 

Grenfell, T. C., and D. K. Perovich, 1984: Spectral albedos of sea ice 
and incident solar irradiance in the southern Beaufort Sea. 
/. Geophys. Res., 89 , 3573-3580, doi:10.1029/JC089iC03p03573. 
Harding, D., P. Dabney, S. Valett, A. Yu, A. Vasilyev, and 
A. Kelly, 2011: Airborne polarimetric, two-color laser altim- 
eter measurements of lake ice cover: A pathfinder for NASA’s 
ICESat-2 spaceflight mission. 2011 IEEE Int. Geoscience Re- 
mote Sensing Symp. (IGARSS), Vancouver, BC, Canada, 
IEEE, 3598-3601, doi:10.1109/IGARSS.2011.6050002. 
Koenig, L., S. Martin, M. Studinger, and J. Sonntag, 2010: 
Polar airborne observations fill gap in satellite data. Eos, 
Trans. Amer. Geophys. Union, 91 , 333-334, doi:10.1029/ 
2010E0380002. 


Krabill, W. B., and Coauthors, 2002: Aircraft laser altimetry 
measurement of elevation changes of the Greenland ice sheet: 
Technique and accuracy assessment. J. Geodyn., 34 , 357-376, 
doi:10.1016/S0264-3707 (02)00040-6. 

Kwok, R., G. F. Cunningham, M. Wensnahan, I. Rigor, H. J. 
Zwally, and D. Yi, 2009: Thinning and volume loss of the 
Arctic Ocean sea ice cover: 2003-2008. J. Geophys. Res., 114 , 
C07005, doi:10.1029/2009JC005312. 

McGill, M. J., D. L. Hlavka, W. D. Hart, V. S. Scott, J. D. Spinhirne, 
and B. Schmid, 2002: Cloud physics lidar: Instrument de- 
scription and initial measurement results. Appl. Opt., 41 , 
3725-3734, doi:10.1364/A0.41.003725. 

, T. Markus, V. S. Scott, and T. Neumann, 2013: The Multiple 

Altimeter Beam Experimental Lidar (MABEL): An airborne 
simulator for the ICESat-2 mission. J. Atmos. Oceanic Tech- 
no l, 30 , 345-352, doi:10.1175/JTECH-D-12-00076.1. 

Perovich, D. K., T. C. Grenfell, B. Light, and P. V. Hobbs, 2002: 
Seasonal evolution of the albedo of multiyear Arctic sea ice. 
/. Geophys. Res., 107 , 8044, doi:10.1029/2000JC000438. 

Schutz, B. E., H. J. Zwally, C. A. Shuman, D. Hancock, and J. P. 
DiMarzio, 2005: Overview of the ICESat mission. Geophys. 
Res. Lett., 32 , L21S01, doi: 10.1 029/2005 GL024009. 

Spinhirne, J. D., 1993: Micro pulse lidar. IEEE Trans. Geosci. 
Remote Sens., 31 , 48-55, doi:10.1109/36.210443. 

Yang, Y., A. Marshak, S. P. Palm, T. Varanai, and W. J. Wiscombe, 
2011: Cloud impact on surface altimetry from a spaceborne 
532-nm micropulse photon-counting lidar: System modeling 
for cloudy and clear atmospheres. IEEE Trans. Geosci. Re- 
mote Sens., 49 , 4910-4919, doi:10.1109/TGRS.2011.2153860. 

Zwally, H. J., and Coauthors, 2011: Greenland ice sheet mass 
balance: Distribution of increased mass loss with climate 
warming; 2003-07 versus 1992-02. /. Glaciol., 57 , 88-102, 
doi:10.3189/002214311795306682. 


